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Coupled CDW and SDW Fluctuations as an Origin of Anomalous 
Properties of Ferromagnetic Superconductor UGe2 
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It is shown that anomalous properties of UGe2 can be understood in a unified way on the 
basis of a single assumption that the superconductivity is mediated by the coupled SDW and 
CDW fluctuations induced by the imperfect nesting of the Fermi surface with majority spins at 
T = 2x(P) deep in the ferromagnetic phase. Excess growth of uniform magnetization is shown 
to develop in the temperature range T < rx(P) as a mode-coupling efi^ect of coupled growth of 
SDW and CDW orderings, which has been observed by two different types of experiments. The 
coupled CDW and SDW fluctuations are shown to be essentially ferromagnetic spin fluctuations 
which induce a spin-triplet p-wave attraction. These fluctuations consist of two modes, spin 
and charge fluctuations with large momentum transfer of the nesting vector. An anomalous 
temperature dependence of the upper critical field Hc2{T) such as crossing of Hc2{T) at P — 11.4 
kbar and P = 13.5 kbar, can be understood by the strong-coupling-superconductivity formalism. 
Temperature dependence of the lattice specific heat including a large shoulder near Tx is also 
explained quite well as an effect of a kind of Kohn anomaly associated with coupled CDW-SDW 
transition. 

KEYWORDS: ferromagnetic superconductor, coupled SDW and CDW fluctuations, UGe2 



§1. INTRODUCTION 

The recent finding of the ferromagnetic superconduc- 
tor UGe2 has been attracting a great interest as the first 
material where both the superconducting (SC) state and 
the ferromagnetic (FM) state are sustained by 5/ elec- 
trons. Contrary to the first expectation, it soon turned 
out that the mechanism of the superconductivity has 
nothing to do with the ferromagnetic spin fluctuations 
associated with FM transition. Indeed, the phase di- 
agram of UGe2 presented in Fig. 1, where Tq and Tgc 
denote the Curie temperature and the SC transition tem- 
perature, respectively, shows that the SC phase is lo- 
cated deep in the FM phase and a steep increase of Tq 
around P = 16.5 kbar suggests the first order transition 
for the ferromagnetic transition.^) This is in marked con- 
trast with the case of the usual unconventional SC phases 
which appear in close proximity to and around the mag- 
netic phase boundary since the superconductivity is me- 
diated by the enhanced magnetic fluctuations around the 
second-order phase transition. The upper critical field 
{Hc2{T)) at T ^ far exceeds the so-called Clogston 
hmit for singlet pairing, Hp{0) ~ 1.84Tsc (Tesla/K),^) 
as seen in Fig. 2. So, it may be natural to confine our- 
selves within a triplet manifold. While phenomenological 
theories based on the ferromagnetism and superconduc- 
tivity appeared recently,^' the mechanism of supercon- 
ductivity consistent with other fundamental properties 
has not been clarified yet. 

The superconducting mechanism appears to be related 
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to another "phase boundary" denoted by T^{P) in Fig. 1 
where a slight decrease of the resistivity has been ob- 
served. As pressure increases, Tx goes down simi- 
larly to Tq and merges into the vicinity of the maximum 
Tgc around P = 11.4 kbar. So, it is natural to consider 
that the fiuctuations related to Tx plays an essential role 
to cause the superconductivity in UGe2 as in the other 
unconventional superconductivity. Although the phase 
boundary was first identified by the small anomaly of 
the temperature dependence of the resistivity p(T), it 
has begun to be recognized quite recently and turned 
out by new measurement that the uniform magnetization 
grows extra (other than the well developed ferromagnetic 
moment) for T < Tx.^'^"^^-' Such metamagnetic behav- 
iors are observed even for P > P^, Px being defined 
by Tx(Px) = 0, under the magnetic field. Indeed, Tx{H) 
reappears under the magnetic field and increases as a rate 
of about 1 K/T near the critical pressure P = P^.^'^^^ 

Furthermore, remarkable anomalies which character- 
ize the superconductivity of UGe2 have been observed in 
Hc2 measurement: As shown in Fig. 2, Hc2 exhibits quite 
unusual T-dependence.^' Namely, iJc2 at T = 11.4 
kbar for which Tsc takes maximum merely shows grad- 
ual increase, while i?c2 at P = 13.5 kbar shows first- 
order like steep increase although Tsc at P = 13.5 kbar 
is smaller than that at P = 11.4 kbar. This behavior can 
be also interpreted as the occurrence of the suppression 
of Hc2{T ^ 0) at P = 11.4 kbar by the magnetic field. 
Such unusual aspects can be understood if the supercon- 
ductivity is induced by the fluctuations associated with 
the transition at Tx as will be discussed in detail in this 
paper. 

So far, the nature of the transition at the phase bound- 



2 



Shinji Watanabe and Kazumasa MiYAKE 



paramagnetism 
^^^^^ 






f e no m.ag neti s m 

T — — ' 

c 












superconductivity 






lujsc fc.^ >^ I 



5 10 15 20 

Pressure [kbar] 



Fig. 1. Temperature-pressure phase diagram of UGe2.'^'' Tc and 
Tgc denote Curie temperature and superconducting transition 
temperature, respectively. The cusp anomaly in the resistivity 
is observed at T = Tx (see text). 



ary T — Tx(P) has not been established yet. However, 
there are some circumstantial evidences that the transi- 
tion may be the coupled SDW and CDW onset due to 
the imperfect nesting of the majority-spin band. 

The first evidence is that the Fermi surfaces of 
the majority-spin band obtained by the band-structure 
calculations is nearly nested both in the a- and c- 
direction.^'^' ^"'^ Although the first stage of calculation 
has been done on the basis of incorrect crystal symme- 
try due to the confusion of experiments, recent calcula- 
tion by H. Yamagami, taking the correct crystal symme- 
try and correct easy-axis (a axis) of magnetization, also 
shows that the Fermi surfaces are mainly composed of 
5/ electrons and have the structure which favors the im- 
perfect nesting in the planes perpendicular to the a- and 
c-axis.-^^^ So, it is expected that the nesting tendency 
is realized in some direction connecting the flat part of 
the Fermi surfaces leading to the enhanced CDW fluc- 
tuations of the majority-spin band around T — T^. It 
is noted that CDW fluctuations of majority-spin band 
is inevitably accompanied by SDW fluctuations. This 
picture is verified by the fact that the ferromagnetic mo- 
ment of the unit cell is almost constant, M ~ l.OfiB, 
along T,. 2. 9. 10) 

The second one is the existence of large lattice entropy 
observed by the specific heat C below T ^ < Tq. In- 
tegrating C/T with respect to T from T = K to Tc, the 
molar entropy is estimated to be 22 J/mole~2.6x A^a^b 
which is far larger than the entropy associated with spin 
degrees of freedom, ln2 x Npjt-Q. The molar entropy due 
to the Fermi quasiparticles is estimated to be 1.7 J/mole 
by means of the Sommerfeld constant limT^o C/T ~ 
0.0314 J/mole-K. This value is much smaller than the 
observed one. The entropy due to the Debye phonon 
also gives only 1/3 of the observed one (see the discus- 
sion in §3). Then the extra entropy should be ascribed 



to the other degrees of freedom. One probable candi- 
date is the softened optical modes coupled to the CDW 
instability of the majority-spin band. 

The third one is the fact that the lattice structure of 
UGe2 is similar to that of a-Uranium in which the CDW 
has been confirmed experimentally by neutron scatter- 
ing in consistent with the nesting vector calculated by 
the band theory. a- Uranium has the zig-zag chain 
arrangement of U atoms along the c axis. UGe2 also 
has the the zig-zag chain arrangement of U atoms along 
the a axis and Ge atoms are located between U atoms. 
Moreover, the T-dependences of p[T) and C(T) of a- 
Uranium are similar to those of UGe2 : The cusp anomaly 
of p{T) and the convex curvature of C (T) at Tcdw have 
been also observed in a-Uranium. Here, one may expect 
that p{T) should raise if the CDW ordering occurs. Of 
course in the case of the perfect nesting, the gap is cre- 
ated at the entire region of the Fermi surface so that the 
system becomes insulator. However, in the case of the 
imperfect nesting, the situation is not so simple. The 
T-dependence of resistivity is determined by the compe- 
tition between two opposite factors: One is the tendency 
that makes the system insulator by the gap creation at 
some parts of the Fermi surface. The other is the effect 
of the reduction of the scattering amplitude among quasi 
particles by the CDW ordering that makes the resistivity 
lower. We remark that NbSe2 where the CDW ordering 
due to the imperfect nesting has been confirmed exper- 
imentally also shows the similar p(T) which decreases 
below TcDW-^''^ 

The fourth one is a rather sharp decrease of the resid- 
ual resistivity p{T 0) across P ^ when P is in- 
creased,^' i*-* implying that the area of the Fermi surface 
recovers as the coupled CDW and SDW ordering van- 
ishes at P = Px- This is because the p{T 0) depend 
directly on the carrier number or the area of the Fermi 
surface and the scattering mechanism is essentially un- 
changed across P = Px, in contrast to the case where T 
changes across with the pressure P fixed. 

In this paper, we show that fundamental properties 
of UGe2 mentioned above can be understood in a uni- 
fied picture that the superconductivity is mediated by 
coupled CDW and SDW fiuctuations due to imperfect 
nesting in the majority-spin band at T = Tx. In §2, we 
show that the growth of extra magnetization for T < Tx 
is understood as a mode-coupling effect of the coupled 
CDW-SDW ordering. In §3, we show that the anoma- 
lously large entropy around T = Tx(P = 0) can be ex- 
plained as contributions from optical phonons softened 
by the Kohn effect due to the onset of CDW ordering 
of the majority-spin band. In §4, we show that the pair- 
ing interaction among quasiparticles in the majority-spin 
band is induced by the ferromagnetic fluctuations as- 
sociated with coupled CDW and SDW fluctuations en- 
hanced around P = P^. Namely, the spin susceptibility 
Xsil ~ 0,iwm) is enhanced by the mode-coupling effect 
of CDW and SDW fluctuations just as these give rise to 
the excess uniform magnetization below T < Tx. These 
fiuctuations induce the spin-triplet "p"-wave paring. It 
is shown that the anomalous temperature dependence of 
the upper critical field Hc2{T) observed at P = 11.4 kbar 
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Fig. 2. Upper critical field for the field applied in parallel to the 
easy axis (a-axis) at various pressures in UGe2.^'^^^ 



and P = 13.5 kbar can be understood by the strong- 
coupling fornialism on the basis of the present mecha- 
nism and the experimental fact that Tx is an increasing 
function of H}'^^ Throughout this paper we use the unit 

such that the Boltzmann constant fee = 1 and the Plank 
constant h = I unless they are explicitly written down. 

§2. Mode-coupling between CDW and SDW 

In this section, we show that the excess growth of uni- 
form magnetization can occur below T < Tx if the transi- 
tion at T = Tx is the coupled CDW-SDW ordering. Such 
a growth of the magnetization in UGe2 has been observed 
by elastic neutron scattering,^) and direct measurements 
of magnetization by ac-method.^'^' -'^"^ In other words, it 
turns out that a central assumption of the present sce- 
nario is supported by a direct experimental evidence. 

An ingredient of such a growth of magnetization below 
T < Tx is the fact that in the free energy there appears a 
mode-coupling term among the uniform magnetization. 
Mo, the order parameter of CDW, Nq, and of SDW, 
M_Q, with a nesting vector Q, in the region T < T^: 

F = Fo + [C{Q,iUn = 0)MoNqM_q + c.c], (1) 

where To denotes the free energy other than the mode- 
coupling term. Here we note that the mode-coupling 
term with uniform charge iVo, Mq and M_q can also ap- 
pear in the free-energy expansion. However, this term is 
considered to be suppressed by the long-range Coulomb 
interaction which works to maintain the charge neutral- 
ity. The mode-coupling coefficient C(Q,ia;„ = 0) in 
eq. (1) can be calculated in terms of the Feynman di- 
agram shown in Fig. 3(a): 



C(Q,iw„ = 0)=5§5^5lQ 

d'P l.cosh-2(£M^). 



(27r)3 4T ^ 2T ^ e(p) - e(p + Q) '^^^ 

where £(p) is the dispersion of the majority-spin band 



and fj, is the chemical potential. The coefficients in 
eq. (2), 5g(> 0), .gQ(> 0), and g'LQ{> 0), denote the cou- 
pling between quasiparticles and spin fluctuations with 
q — 0, charge fluctuations with q = Q, and spin fluctu- 
ations with q = — Q, respectively. 

First, we consider the perfect-nesting case. The per- 
fect nesting occurs in the bipartite lattices at half filling 
with /i = 0. In this case, the relation e(p -|- Q) = — e(p) 
holds and the mode-coupling coefficient is given as fol- 
lows: 

^2T> 



C(Q,ia;„ = 0)= ^offqff-qf^ / ^^^(0^^ cosh" 



- 9o9q9-q 



1^=0 



(3) 



where D{^)9. = J (fp/{2wfS{e{p)-iJ,-0 is the density 
of states measured from at the Fermi level, fi, with Q, 
being the volume of the system. The last approximate 
equality in eq. (3) is valid when Tx is assumed to be much 
smaller than the bandwidth of quasiparticles. Therefore, 
the coefficient C vanishes as far as the D{^) is symmetric 
function of ^ at ^ ~ 0. This result can be understood 
as follows: As is well known, the number density of the 
majority spin band, under the fixed chemical potential, 
remains unchanged even under the growth of the coupled 
CDW and SDW ordering, Ac-s cx MqNqM^q, in such 
a case. This can easily be seen if one considers the one- 
dimensional simple lattice as an example in Fig. 3(b): 
By the CDW-SDW transition the rearrangement of the 
charges with majority spins occurs but the net magneti- 
zation does not change as shown in Fig. 3(c). 

In case that D{(_) has a weak asymmetric part with 
respect to ^ = 0, the coefficient C does not vanish so 
that the number density of majority spin band changes 
in general. However, it may be reasonable to assume 
T'(^) is symmetric around ^ = in the case where the 
perfect nesting is realized as in two-dimensional(2D) and 
three dimensional(3D) Hubbard model with only near- 
est neighbor hopping. In one-dimensional case, how- 
ever, D{(_) has a weak asymmetry in general under the 
condition of perfect nesting. In such a case, the sign 
of the magnetization change cannot be determined by 
the sign of the term of 0(AQ_g), but is determined 
by the sign of the term with logarithmic correction, 
AQ_glnAc~s- In any case, the effect is small one pro- 
portional to \dD{S^)/d^\^=o\ <C -D(^ 0), in contrast to 
the case of inperfect nesting where the effect is propor- 
tional to D(^ = 0) as will be discussed below. 

Next we consider the imperfect-nesting case. In this 
case, the energy gap opens at the Fermi surface on 
the zone boundary transferred by the nesting vector Q. 
Hence, the coefficient eq. (2) does not diverge but has 
a finite value. The important point is that the sign of 
C{Q,iu}n = 0) changes according to the relative posi- 
tions between the Fermi surface and the nesting zone 
boundary. 

As an example let us consider the majority-spin band 
on the two-dimensional square lattice 

e(p) = -2t{cos{px) + cos(py)} - 4t'cos(px) cos(py), 
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with the nearest-neighbor hopping t and the ncxt- 
nearest-neighbor hopping t' . In Fig. 3(d) and Fig. 3(c) 
the Fermi surface with fi = and nesting zone boundary 
is plotted by the sohd hne and the dashed hue, respec- 
tively. 

For positive t' the Fermi surface with /i = is located 
inside the nesting zone boundary as in Fig. 3(d). In this 
case the mode-coupling coefficient has a negative value 
C{Q,iuJn = 0) < 0. Namely, in this case the uniform 
magnetization Mq increases for T < T^. This can be 
understood as follows: The energy gap opens at the part 
of the Fermi surface along the nesting zone boundary, 
which is drawn by the shaded spots in Fig. 3(d). By the 
gap creation, the energy band inside the zone bound- 
ary shifts downward. As for the area at which the gap 
opens the situation is the same as the perfect nesting 
case. Namely, the shaded spots in Fig. 3(d) do not con- 
tribute to the change of uniform magnetization. How- 
ever, the important difference from the perfect-nesting 
case is that the Fermi surface which deviates from the 
nesting zone boundary (i.e., the part where the gap does 
not open) shifts downward by the gap creation. Hence, 
the electrons flow into the majority-spin band from the 
minority-spin band so that the uniform magnetization 
increases. 

On the contrary, for negative t' the Fermi surface with 

= is located outside the nesting zone boundary as in 
Fig. 3(e). In this case the energy band shifts upward by 
the gap creation and it is shown that the mode-coupling 
coefficient has a positive value C(Q,iLi;„ = 0) > 0. 
Namely, in this case electrons flow into the minority-spin 
band from the majority-spin band so that the uniform 
magnetization Mq decreases for T < Tx. 

Here we have discussed the change of the uniform mag- 
netization referring to the one- and two-dimensional lat- 
tices as typical examples to draw a general conclusion. 
In the real systems, the shape of the Fermi surface is 
often more complicated, but our statement derived here 
can be applied to the Fermi surface transferred by the 
nesting vector Q. As for UGe2, the Fermi surface of 
the majority-spin band is located inside the nesting zone 
boundary according to the recent band-structure calcu- 
lation. Namely, it seems to correspond to the case 
(d) in above examples, which should exhibit the increase 
of the uniform magnetization for T < Tx as observed 
experimentally .2.9. lo) 

Concluding this chapter, let us remark that the mech- 
anism of magnetization changes due to simultaneous or- 
dering of CDW and SDW with the equal nesting wave 
vector is quite general one which may be found in other 
systems with such a property. It is also remarked that 
the sign of magnetization change depends on the cur- 
vature of the Fermi surface near the nesting vector as 
discussed above. 



§3. Specific heat due to Kolin anomaly 

In this section, we discuss how the temperature de- 
pendence of the speciflc heat C(T) at T < Tc, especially 
a gentle hump around T = Tx, is understood from the 
present point of view that Tx is the onset temperature 



C(Q, m=0) =0 




-It -71/2 Jt/2 71 
C(Q, m=0) < 



-TC -71/2 71/2 7C 

CiQ, /CO, =0) > 



Fig. 3. (a) The mode-coupling term between the CDW and the 
SDW in the lowest order. The band picture (b) and the real- 
space picture (c) by the perfect nesting in the majority-spin band 
in the one-dimensional simple lattice. The Fermi surface with 
fi = (solid line) and the nesting zone boundary (dashed line) 
in the two-dimensional square lattice in the majority-spin band 
e(p) = -2{cos(p^) + cos{py)} — it' cos{px) cos{py) with (d) t' = 
0.4 and (e) t' = —0.4. The imperfect nesting by vector Q = 
(tt, it) causes the gap at the shaded spots. 



of coupled CDW-SDW ordering. The specific heat due 
to the Debye phonon cannot afford to explain the global 
behavior of C(T)/T as can be seen in Fig. 4 in which the 
Debye contribution Cd/T is compared with the experi- 
ment at ambient pressure. The Debye contribution is 
given as 



^ = A 
T 



4 

2^2 



7^:77 / duj- 



1 



with the coefficient A defined by 



A = 



9Nk 



(4) 



(5) 



D 



where N and ujy> denote the unit-cell number included in 
the system and the averaged Debye energy, respectively. 
In the low temperature limit, T ^ wd, eq. (4) is reduced 
to the conventional form as 



T 



5^3 



(6) 



Comparing the experimental value for the coefficient of 
T^ term of eq. (6), 3.2x10-* J/K'^-mole,^) we obtain the 
Debye energy as wd = 182 K. Using this value for wd, we 
have ploted Cd/T in Fig. 4. In order to account for the 
hump around Tx, there should be an extra contribution 
to the entropy from other than the magnetic origin be- 
cause the difference between Cd /T and the experimental 
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value is far larger than the jump of C/T at T = Tc, the 
ferromagnetic transition. 
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Fig. 4. C/T vs T. The solid circles represent the experimental 
values,^' the solid line the theoretical ones, which consist of the 
Debye contribution (dotted line), the cotributions from optical 
modes with the Kohn anomaly (dashed line), and the fermionic 
contributions from quasiparticles (solid line). 



The optical phonons are subject to the Kohn anomaly 
if the Fermi surface of quasipartices is nested as expected 
in the majority-spin band of UGe2. We parameterize 
the dispersion of optical phonons exhibiting the Kohn 
anomaly as follows: 



e(k) = ^[s\\k.f + sl{kl + klf + . 



(7) 



where the wave vector k is measured from the nesting 
vector Q, and we here take the direction along which 
the nesting occurs as fc^-direction without loss of gener- 
ality. The gap A parameterizes the degree of softening 
of the phonons, and s\\ and s± denote the sound veloc- 
ities for longitudinal and transverse directions, respec- 
tively, with respect to Q. There exist two independent 
optical-phonon modes corresponding to two Ge ions in 
the unit cell. For simplicity, we assume that the disper- 
sions of two modes are identical. Then, the energy due 
to optical-phonon modes with Kohn anomaly is given as 



£(cJ2,Lj) 



1' 



(8) 



The 



where Wz = s\\kz, w = s±^kl + k^, S = s±/s\\. 
Debye energy lo^ is different from cjd for the acoustic 
modes in general. 

In the high temperature limit, T ^ lu^, eq. 
reduced to 



(8) is 



1 



(9) 



Jo Jo ^ 

The coefficient S is determined so as to reproduce the 
value of Dulong-Pctit: 

„ 6A^fcB 



UJ- 



*3 



X 2. 



(10) 



The specific heat Ck is obtained from eq. (8) as follows: 



Ck 
T 

1 

'2T~dT 



1 

dw. 



dujz 



dujuj 



1)2 
1 



c (e^/^-l)2 J-(^^) 

The second term crucially depends on the temperature 
dependence of the gap A defined as the lowest optical- 
phonon energy in eq. (7) . To discuss the relationship be- 
tween the temperature dependence of the specific heat 
and the effect of the electron-phonon coupling, we show 
the Dyson equation for the phonon Green's function 
D{c[,iujn) in Fig- 5(a). Here Do{q,iujn) denotes the 
bare Green's function for phonons. 11 and Smc denote 
the electron polarization function and the selfenergy for 
the mode coupling, respectively. The softening of the 
optical-phonon mode with gap A at the wave vector Q 
which should be the same as the nesting vector is illus- 
trated in Fig. 5(b). On the mean- field level where Smc 
is set to be zero, the T-dependence of A is given near 
the ordering temperature like 



J\T - T, 



MFi 



(12) 



where J is a positive constant which is determined by 
the band structure of quasiparticles and the electron- 
phonon coupling constant. Then, Ck/T given by eq. (11) 
exhibits a discontinuous jump at T = through the 
factor dA^/dT as shown in Fig. 5(c). However, in case 
the mode-coupling term cannot be neglected, in the limit 
A ^ the mode-coupling effect gives the self-energy of 
the relevant optical phonon mode the term like —Mi + 
M2A, Ml and M2 being positive constants with weak 
temperature dependence near as discussed in ref.^^^ 
(see Fig. 5(d)). Thus the self-consistent equation for A 
takes the following form at T > Tx, Tx being the ordering 
temperature determined self-consistently as below: 



J(T - T, 



MF\ 



Ml - M2A. 



(13) 



By solving this equation near 
obtain 



MP 



(Mi/J), 



we 



J'^{T ~Ty_f /Ml 



(14) 



For T < Tx, similar temperature dependence is obtained 
by solving much more complicated equation due to the 
existence of ordering. Namely, it turns out that in case 
the effect of the mode coupling is relevant, the jump 
in the specific heat at the transition temperature Tx is 
smeared out by the fluctuations. We consider that UGe2 
may correspond to such a case where the mode-coupling 
term cannot be neglected for the fluctuations due to the 
strong electron-phonon coupling. Thus, to calculate the 
specific heat of UGe2 we introduce a parameterization 
for T-dependence of A as follows: 

This is reduced in low- and high-temperature limits as 
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follows: 



A2(T) 



/\2 

a 



(16) 



D D„ O n ^ A' 

I [ Mode Coupling 

D{q = D^lq M,) " n(<? ,icq,) - £mc 



= aliq) - \n(q,0)\ + Emc 




Mode Coupling 




Jump in C at 7;; No jump in C at 7i 



Fig. 5. (a) Dyson equation of phonon Green's function including 
the selfenergy of the mode coupling, (b) Softening of optical 
phonon mode at the nesting vector Q. Temperature dependence 
of the gap in the mean-field framework (c) and in the framework 
including the mode-coupling term (d). 



In Fig. 4, the contribution from eq. (11) is shown as 
dashed line for parameters lo^=21Q K, Ao=70 K, a=0.5, 
and (5=1.0. One solid line in Fig. 4 is a contribution from 
the Fermi quasiparticles 7quasi which is determined ex- 
perimentally as 7quasi = limy^o C/T ~ 0.0314 J/mole 
K^.^) The other solid line is for Cth_/T which is the sum 
of C-ajT, Cy^IT, and 7quasi- The curve Cth/T repro- 
duces quite well the global behavior of observed C/T. In 
Fig. 6, C/T (solid circle), Cth/ (solid hue), and their dif- 
ference AC/T = C/T — Cth/T (solid square) are shown. 
The origin of /S.C /T near T = Tc is apparently the mag- 
netic entropy of local component of magnetization, while 
that of a hump around T=30 K may be attributed to 
the mass enhancement due to the criticality of coupled 
CDW-SDW ordering at T = ~35 K, where the en- 
hancement of a few times of 7quasi is expected as observed 
at (T,P)~ (0,P,).6) 

If we had not taken into account the effect of Kohn 
anomaly, we could not account for the gentle hump of 
C jT around r=35 K by the Einstein model which simu- 
lates the two optical phonons with frequency we. Indeed, 
its contribution Ce/T is given as 

6A^feB WEe'^^/^ ^^^^ 



Ge 
T 



T3 (go^E/T _ 1)2 



This has a rather sharp maximum C^'^^/T ~ 76/ci;| J 
mol~"'^K~^ at T ~ O.SOt^E; so that it is difficult to re- 
produce a broad hump structure of C /T observed by 
experiments. 




Fig. 6. Difference of C/T between experiment and model. The 
solid circles represent the experimental values,^) solid line theo- 
retical ones of model, and the solid squares the differences. 



§4. Pairing Interaction and Strong-Coupling Ef- 
fect on Hc2 

As we discussed in § 2, the excess ferromagnetic com- 
ponent begins to arise at T = Tx, where the uniform 
spin susceptibility is diverging due to the mode-coupling 
effect among charge and spin fluctuations with the nest- 
ing vector Q. While the charge susceptibility can be 
also diverging at T = due to the mode-coupling ef- 
fect among the uniform charge A^q=o and two modes of 
SDW, Mq and M_q, uniform charge fluctuations are 
suppressed by the long-range Coulomb interaction which 
works to maintain the charge neutrality as mentioned in 
§ 2. Then, the pairing interaction between quasiparti- 
cles with majority spin, due to those mode-coupling ef- 
fects, is mediated mainly by the coupled two modes of 
the spin susceptibility Xs and the charge susceptibility Xc 
as shown in Fig. 7(a), but not by the coupled modes of 
Xs and Xs- The coupled fluctuations of Xs and Xc which 
are related to the excess ferromagnetic moment are re- 
garded as the ferromagnetic spin fluctuation Xs which is 
enhanced for the small momentum transfer, i.e., q '--^ 0. 
The important point is that the coupled two modes can 
make the scattering wave vectors for the Cooper chan- 
nel, q = k — k', zero which is profitable for the triplet 
pairing (see Fig. 7(a)), in marked contrast to the single- 
mode case where Xs or Xc with finite Q gives a small 
repulsive part for the triplet pairing as has been well 
recognized. ^""^'^^ 

These coupled fluctuation propagators give non-zero 
contributions even for the external wave vector g = and 
the bosonic matsubara frequency 7^ 0, '^m = ^rmrT. 
If we parameterize these propagators as 



Xs(Q + q,iwm) ~ Xc(Q + q, it^™) 

Xo 



(18) 



770 + Aoq"^ + Co\uj^ 
the propagator Xs of the coupled fluctuations is given as 



Xs (g, it^m ) = r 2_, Xs (Q + q' , 
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XXs(Q + q',t^m') 



q' 



dx . X 
— cotn — 
TT 2T 



q+q' 



■ Co\uJm\) 



■ ip^, + X^ {(fq+q' 



Ca\LUrn\)^ + X^ 



where ipq = 770 + Aqq^. It is easy to verify 
that Imiq^Q Xs{q,iLL!m) is finite (non-zero) and con- 
tinuous as a function of a;,„ and diverges in pro- 
portion to Txoi4:nAo)~^ /rjo in the 2D system, and 

^Xo(8'''^o^^)~^/\/^ ill the 3D system, as approaching 
the onset of coupled CDW-SDW ordering at finite tem- 
perature, i.e., LOm =^ and 770 — > 0. Near the quantum 
critical point at T — 0, these divergent behaviors are 
given as ln(l/r7o) in the 2D system, and (const. — -^770) 
in the 3D system, respectively. It is noted that the latter 
expression in the 3D system exhibits a sharp peak with 
cusp at ?7o, while it does not diverge there. 

By careful inspection of the expression of right hand 
side of eq. (19), one can see that the main depen- 
dence of Xs(9:i^m) on q and cj,„ arises through a 
combination Aqq^ + Co\uJm\- An explicit dependence 
of Xs{q,iuJm) on Agq'^ + Co|wm| is not simple com- 
pared to those for Xs and Xc eq. (18), and depends 
on the space dimensionality. Indeed, the main part 
of Xs{qA^m) is proportional to TxIC^ttAq)^^ {Aoq-^ + 
Co|w,„|)-i ln[(2r7o + Aoq^ + Cq\u;^ \)/2t]o] in the 2D sys- 

tem, andTxliiTTAl/Y'/iVm + Vm + i^oq^ + Co|t^™|)/2] 
in the 3D system. The corresponding expressions at 
T = near the quantum critical point are too compli- 
cated to be written here. 

Since the purpose of the present paper is to give a 
theoretical basis for qualitative but unified understand- 
ing of the superconducting mechanism of UGe2, espe- 
cially an origin of anomalous temperature dependence of 
Hc2{T), we adopt the following phenomenological form 
for Xs{q,iuj„i) as 

Xo 



there remains an incoherent contribution which arises 
from high energy processes and cannot be described by 
the response function of quasiparticles. It is quite cru- 
cial for the behavior of incoherent processes whether the 
relevant quantity is conserved or not, while it is not so se- 

(19) 

vere for the quasiparticle response which mainly reflects 
the structure of the phase space of low energy particle- 
hole continuum of quasiparticles as far as the effect of 
the impurity scattering is neglected. 

On the basis of the reasoning above, the pairing in- 
teraction V{q,iuJm) for q = k — k' ~ is postulated 
as 



y(q,iw™) ~ f7- |rpxs(q,iw™), 



(21) 



where U is the short-range (q-independent) repulsive in- 
teraction, I* is the coupling between quasiparticles and 
the coupled CDW-SDW fluctuations. The wavenumber 
dependence of / is assumed to be neglected around g '--^ 0. 
The pairing interaction arising from the exchange of the 
single mode of SDW and CDW fluctuations with the 
nesting vector Q gives a small repulsive part for the 
triplet p-wave interaction as mentioned above. However, 
this contribution is not considered to play a crucial role 
because the wavenumber dependence of Xs(q) and Xc(q) 
is weak at q ~ 0. So, we neglect its contribution for the 
moment. 



-k', -;£,', k-k'~ 

%^ik-k\ ie>-i£^;) 

III 

2 Xo 




1;] + A(k-ky+ C {£,,- 



(b) 



-k, -ie„ 



Xs(q,iwm) 



IT + Aq^ + C\uj^\' 



(20) 



This reproduces approximately a characteristic behavior 
of Xs('3', iwm) near the critical point of coupled CDW- 
SDW ordering if the parameters 77, A, and C are con- 
stants different from 770, Aq and Co in eq. (18) and deter- 
mined so as to reproduce key experimental observations. 

The q-dependence of the vertex connecting the cou- 
pled CDW-SDW fluctuations and the quasiparticles de- 
pends on a dynamical structure of the fluctuations. If the 
magnetization is the conserved quantity, the vertex for 
q ~ and ujrn 7^ should vanish as can be verified by a 
perturbational calculation of the triangle diagram shown 
in Fig. 7(a). However, if the magnetization is not con- 
served due to the strong spin-orbit interaction at U-sites 
as in the present case, the vertex does not vanish even 
for q = and Um ^ 0. The mode mediating the pairing 
interaction shown in Fig. 7(a) is that arising from an in- 
coherent process in the sense of the Fermi liquid theory. 
It is known that such a contribution vanishes in general 
if the relevant quantity of the fluctuations is conserved 
one.^'*-' On the other hand, if it is not conserved one, 



Fig. 7. (a) Paring mechanism mediated by coupled CDW (xc) 
and SDW (xs) fluctuations with the nesting vector Q. (b) Self- 
energy of electrons due to one-fluctuation exchange process. 



For simplicity of calculations, we assume a three- 
dimensional spherical band with majority spins. Then 
the triplet p-wave component V{£ — 1) of V^(q, iwm) for 
q = k — k' with |k| = |k'| = k-p is given by 

y(£= l,iw™) « -irp / — sin0V3cos0Xs(g,i'^m), 
Jo ^ 

(22) 

where q = 2A:f| sin 6*72 1, with cos6' = (k • k')/A:|. Sub- 
stituting eq. (20) for Xs and after simple algebra, we 
obtain 



V^(^ = l,ic^,„) « - 



V3|rpxo 
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-1 + 



B 



In 



{B + 2) 



B 



where w-dependent B is defined as 



(23) 



(24) 



The corresponding selfenergy S of electrons due to the 
one- fluctuation exchange process (see Fig. 7(b)) is given 
as follows: 



I](p,ie„)=T^^|r 

e„/ p' 



*i2Xs(p - P',ien -ien' 



ie„ 



, (25) 



where e„ = (2n + 1)-kT is the fermionic Matsubara fre- 
quency. The summation in cq. (25) with respect to p' 
can be separated approximately and reduced to 

E(p,ie„)=T5:|rr^^^ 

/■TT 

J — sin 6'Xs (g, ie»i - ie^' ) , (26) 



X 



where q = 2fcF|sin^/2|. Here we have approximated as 
|p| ~ Ip'I ~ fcp in Xs in eq. (26). Then, the renormalized 
frequency e„ = e„ + iS(ie„) is given approximately as 
follows: 

|/*|2^oiVFsign(e„0 



En = + 7I"T ^ ■ 



X In 



B(e« -£„') + 2 



i?(e„-6„0 ' ^^^^ 

where N-p is the density of states of band electrons. In de- 
riving eq. (27), we have discarded the real part of S(ie„) 
on the understanding that it can be absorbed into the 
chemical-potential shift. 

Here we introduce the function for renormalization ef- 
fect as 

r + + 2^;* 



In 



M(uJm) = A- 



r+ \u)n 



In 



T + 2E* 



where A, F, and E* are defined as 



T + 2E* 



and 



C" 



C 



(28) 

(29) 
(30) 

(31) 



Here F is the damping rate of coupled CDW and SDW 
fluctuation and E* is the order of the unrenormalized 
Fermi energy of band electrons. 
Then, eq. (27) is written as 

CX) 

e„ = e„ + ttT ^ M(e„ - e„)sign(e„). (32) 

m— — oo 

After some algebra, we can show that the left hand side 



of eq. (32) is reduced, in the limit of T — > and up to 
the linear term of e„, to 

e„ = e„(l + A). (33) 

Namely, A is the coupling constant giving the mass en- 
hancement factor. Of course, there may exist other con- 
tributions to the mass enhancement arising from single- 
mode processes of exchanging spin and charge fluctua- 
tions with large nesting vector Q. However, this contri- 
bution is expected to give only moderate enhancement 
for the Sommcrfcld constant if the fluctuations are dom- 
inated by the nesting as has been discussed elsewhere. 
So, we neglect this contribution for simplicity of analysis. 

In terms of A. F and E* , the pairing interaction of 
p-wave (€=1) is written as 



V3 



2A 



In 



T + 2E* 



-1 + 



E* 



2E* 



In 



2E* 



(34) 



On the derivation above, the pairing interaction V for 
a relevant channel can be expressed in a separable form 
in the wave vector and frequency space as 

V(k - k',ie„ - ie„0 = 0(k).^(k')w(ie„ - ie„0, (35) 

where (/'(k) is the wavefunction of relevant pairing chan- 
nel, and V represents the frequency dependence of inter- 
action due to the retardation. In the present case, 0(k) ~ 
VSkz, and v is given by eq. (34). Thus, the gap function 
A also has the separable form A(k, ie„) = (/)(k)A(ie„) 
and satisfies the gap equation 



A(ie„) = -r^^t;(ie, 



') 



|0(k')|2A(ie„O 



"^^' + e^ + |0(k')A(ie„OP"' ^^^^ 

where is the dispersion of quasiparticles without the 
renormalization due to the many-body effect of eq. (26), 
and e„ is the renormalized frequency due to the many- 
body effect of cq. (27). Explicit form of e„ in the present 
model is given by eq. (32). 

The linearized gap equation is then given as follows: 

A(ie„) = -TY^vi^n - i6„0 A(ie„0 E 



(37) 



At this stage, it may not be so bad to approximate 
|(/ii(k)p as an average (|(/)(k)p) near the Fermi level. This 
quantity is nearly equal to unity due to the normaliza- 
tion condition. Namely, eq. (37) can be approximated 
by the same one for the conventional s-wave pairing. 

So far, we have constructed the general framework 
of the pairing mechanism mediated by coupled CDW 
and SDW fluctuations. To discuss the real materials 
rigorously, it is necessary to take into account of the 
band structure and the Fermi surface, accurately. Here- 
after, we would like to discuss the anomalous behaviors 
of Hc2{T) under pressures in UGe2. However, there is 
not enough information about the band structure and 
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the Fermi surface responsible for the superconductivity 

under pressure and magnetic field in UGe2 at present. 
Hence, we proceed the discussion on the basis of the 
above framework with three-dimensional spherical band 
and consider the features of UGc2 by taking into ac- 
count of the most-important factors such as pressure and 
magnetic-field dependences of mass-enhancement factor 
and Tx as will be discussed below. It is an important task 
to consider the accurate band structure and to determine 
the pairing symmetry of the orbital part in the next stage 
of the theoretical approach. However, we think that the 
anomalous behaviors of Hc2 under pressures do not de- 
pend on such details but are caused by the key factors 
which will be explained below. 

Let us summarize again the behaviors of upper-critical 
field under pressures (see Fig. 2) and give its intuitive un- 
derstanding. At P — 11.4 kbar, iJc2 shows only gradual 
increase in spite that Tgc was originally high at -ff = 0. 
This is because Tx, around which the attractive interac- 
tion is induced, goes up and the enhanced fluctuations 
associated with Tx fades away as H increases. On the 
other hand, at P = 13.5 kbar Tsc is less than the maxi- 
mum value at i? = 0. However, as H increases, Tx goes 
up to be close and the attractive interaction is enhanced 
so that Hc2 grows rapidly. Then, Hc2 at P = 13.5 kbar 
crosses above that of P = 11.4 kbar. As for P = 15.3 
kbar it is interpreted that i?c2 reaches T = before Tx 
closes up under magnetic field and thereby shows only 
the gradual increase. 

The strong-coupling formalism for Hc2 originally de- 
veloped for s-wave superconductors on the; c;k;c;tron- 
phonon mechanism, ^^'^^^ can be extended rather easily 
for the present case. It is because the effect of magnetic 
field can be incorporated into the theory just as in the 
conventional case as far as the vector potential only cou- 
ples with the center of mass coordinate of the wave func- 
tion of the Cooper pair, but does not affect its relative 
coordinate. This is the case in usual situations since the 
size of the Cooper pair is far smaller than the width of 
the Landau orbital except for H ~ -ffc2(T = 0). Thus we 
can use the strong-coupling formalism for i/c2 developed 
for treating the s-wave pairing on a slight modification. 

Starting from the Eliashberg equations, eq. (32) and 
eq. (36), one obtains the linearized gap equation under 
the magnetic field H as follows: 

A(«en) = ^Y^ - em)J^(em)A(ie„), (38) 

m 

where the effective coupling V is given by 

V{um) = -N-pV{e = 1, iWm), (39) 

where V{£ = 1, iUm) is defined by eq. (34). The integra- 
tion kernel K in eq. (38) is given as 

where e is the elementary charge, vp is the Fermi veloc- 
ity of band electrons and e„ is defined by eq. (32). In 



the present case where only the electrons with majority 
spin participate in the formation of the SC state, the 
non-unitary triplet pairing is expected. Then the para- 
magnetic effect due to the Zeeman splitting has been 
neglected in deriving eq. (40). The strong-coupling ef- 
fect is included in the frequency dependence of the gap 
function A(ie„) in eq. (38) and e„ in eq. (40). The upper 
critical field Hc2 is determined by the condition that the 
linearized gap equation starts to have a nontrivial solu- 
tion A(ie„) 7^ as iJ is decreased with T fixed. The 
resultant equation to be solved is 



det 



T ~ 

Snm - irr^i'^n - em)-^(em) 
iVp 



= 0, (41) 



where 6nm is the Kronecker delta symbol. 

Our strategy for discussing Hc2{T) is as follows: First, 
we determine the pressure dependence of A in eq. (28), 
giving the mass enhancement factor, by experimental 
data of specific heat.^^ Second, we estimate the pres- 
sure dependence of parameters F, eq. (30), and E*, 
eq. (31), so as to reproduce Tsc- Third, we deduce the 
magnetic-field dependence of the parameters with the 
use of a functional form of T^{P,H) experimentally ob- 
tained.9'1^) Then, H^2{T) is calculated from eq. (41). 

First of all, we discuss the P-dependence of A at 
H = 0. In Fig. 8 open circles with an error bar de- 
note the ratio of the specific-heat coefficient under finite 
and zero pressures, 7(P)/7(0), which are measured at 
H = 0.^-' This ratio can be written in terms of the mass- 
enhancement factor as7(P)/7(0) = [l + X{P)]/[l + X{0)]. 
The enhancement of 7(P) for P > 13 kbar is considered 
to be contributed by two factors: one is the fluctua- 
tion associated with the anomaly at Tx and the other 
is the critical FM fluctuation near Tq which we do 
not consider here as irrelevant for the superconductiv- 
ity. We consider the former is essential for inducing the 
SC state. Hence, we introduce a biquadratic function of 
[1 + A(P)]/[1 + A(0)] around P = Pc ~ 13 kbar to fit 
7(P)/7(0) as the dashed line in Fig. 8, although experi- 
mental data for 7(P) remains to increase gradually even 
in the P > 13 kbar region. We assume that the mass en- 
hancement at zero pressure gives no contribution to the 
SC state for 10 < P < 17 kbar, i.e., we set A(0) = 0. 
In this way, we now obtain the pressure dependence of 
the mass- enhancement factor, A(P) in eq. (28), which is 
denoted by the shaded line in Fig. 8. 

Secondly, we determine the P dependence of F, 
eq. (30), and E* , eq. (31), so as to reproduce Tsc{P) 
at = by using A(P) determined as above. In Fig. 8 
open triangles denote 100 times of the SC transition tem- 
perature, i.e., lO^Tsc- We assume the shape of r(P) as 
the dash-dotted line in Fig. 8, which becomes narrow 
toward the critical pressure P ~ 13.5 kbar. The open 
squares denote the values of E* which reproduce the ex- 
perimental data of Tsc at each pressure. Here, we intro- 
duce a biquadratic function of E* (P) as the solid line in 
Fig. 8 to fit discrete data denoted by open squares. Then 
we obtain the P-dependence of A, eq. (29), F, eq. (30), 
and E* , eq. (31) from the experimental data of 7 and 
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Fig. 9. (a) Expcriiiiciital data of Tx as a function of filed applied 
parallel to the crystallographic a-direction at P = 11.4 kbar 
and P = 13.5 kljar.^*^' (b) Introduced function of T^{P,H) to 
calculate H^2 at each pressure. 



Fig. 8. Pressure dependence of parameters at H = G. Open cir- 
cles with an error bar denote the ratio 7(P)/7(0) measured with- 
out magnetic field. Open triangles denote ten times of the 
SC transition temperature Tgc-'^^ Open squares denote E*{P) 
which reproduce Tgc by using A(P) and r(P). The dashed line 
and the shaded Une denote [1 + A(P)]/[1 + A(0)] and A(P), re- 
spectively. The dash-dot line and the solid line denote r(P) and 
E*, respectively. 



Tsc in the H = case. If we choose another shape 
of r(P), the shape of E*{P) also changes to reproduce 
Tsc- However, the spectrum of coupled CDW and SDW 
fluctuations which mediate the superconductivity is con- 
sidered to be sharp at Tx and thereby we here adopt the 
function T{P) as shown in Fig. 8. 

Thirdly, we discuss the i7-dcpendcncc of the param- 
eters A, r, and E* . Here we note the if-dependence of 
Tx is identified by the rapid change of the slope of the 
resistivity p(T) under the magnetic field at each pres- 
gyj.g9,28) jg shown in Fig. 9(a). We see that Tx 
depends on H linearly at P = 11.4 kbar and increases 
steeply at P = 13.5 kbar. Incorporating these behaviors, 
we introduce a function T^{P,H) as shown in Fig. 9(b). 
The key point is that Tx is an increasing function of H 
and a decreasing function of P. By using this form of 
function T^{P,H), the H dependence of A, eq. (29), and 
r, eq. (30), and E* , eq. (31), is taken into account in the 
following way: We can obtain the pressure P at iJ = 
which gives the same value of Tx at P and H by using 
the relation 

T4P,H = 0)=T^{P,H). (42) 

This always gives P < P for > 0. Namely, applying 
the positive magnetic-field turns out to be equivalent to 
the negative pressure at H = 0. In this way we can 
take into account the effect of the magnetic field on pa- 
rameters A, F and E* which appear in the gap equation 
eq. (41) of our model. 

On these estimation of parameters, we can calculate 
Hc2 at each pressure as follows: The pressure P and 
the magnetic field H are fixed. With use of eq. (42), 
the reduced pressure P = P is determined. Then, the 
linearized gap equation eq. (41), with parameters A(P), 
r(P) and E*{P) introduced in Fig. 8, is solved to obtain 
Hc2- The Fermi energy ep and the Fermi velocity vp of 



unrenormalized band electrons have been set to repro- 
duce the order of magnitude of observed Hc2 (T ^ 0) . 

Results of iJc2(T) for P = 11.4, 13.5 and 15.3 kbars are 
shown in Fig. 10. At P = 11.4 kbar, where TsciH = 0) is 
nearly maximum, Hc2{T) merely increases gradually and 
monotonously as T decreases. This is because applying 
the magnetic field gives negative pressure which makes 
Tsc{H) decreased, as seen in Fig. 8. This reflects the 
fact that Tx(iy) goes up and away from Tsc{H = 0) as 
H increases at P = 11.4 kbar. On the other hand, at 
P = 13.5 kbar, where Tsc{H = 0) is rather smaller than 
the maximum at P = 11.4 kbar, Hc2{T) increases steeply 
and shows first-order like "transition" at = 1.8 T. This 
anomalous increase is understood as a result of sudden 
appearance of Tx at i7 = 1.8 T due to the magnetic field 
as shown in Fig. 9(a). Alternatively, such an increase can 
be understood by the fact that the applied magnetic field 
has an effect of negative pressure which makes Tsc{H = 
0) shift to the maximum value at P = 11.4 kbar from 
P = 13.5 kbar as shown in Fig. 8. From these results 
we can see that the crossing of i?c2(T) at P = 11.4 kbar 
and P = 13.5 kbar is caused by two counter effects of 
the magnetic field: The suppression of Hc2{H) at P = 
11.4 kbar and the enhancement of Hc2{H) at P = 13.5 
kbar. At P = 15.3 kbar Hc2 shows the gradual and 
monotonous increase as T decreases. It is because the 
magnetic field H < Hc2{T = 0) is not enough to recover 
Tx{H) > so that the concave T-dependence of ifc2(T) 
is not realized. 

If the H dependence of Tx at P = 13.5 kbar is grad- 
ual rather than the square-root-like behavior used in the 
present calculation as in Fig. 9(b) (for example, if Tx has 
a linear dependence with respect to H), the first-order- 
like jump in Hc2 does not appear. However, i?c2(T) 
still increases steeply and crosses that of P = 11.4 kbar. 
Then, the detail of the shape of i?c2(T) changes accord- 
ing to the H dependence of Tx, but the global feature of 
Hc2{T) under pressure such as crossing does not change. 

The above result of i?c2(T) is qualitatively the same 
as that obtained on a more simplified phenomenological 

model in which the effective coupling X{uj,m), eq. (29), 
and the function for renormalization effect Miojm.), 
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Fig. 10. The upper critical field _ffc2 calculated based on the 
strong-coupling-superconductivity formalism explained in the 
text at P = 11.4 kbar, 13.5 kbar and 15.3 kbar. 



eq. (28), are approximated as 

A(c^™)«MKO« (43) 

As shown in ref,^^'^®-* this model explains the crossing of 
Hc2{T) at P = 11.4 and 13.5 kbars when the parameters 
A and F are fixed in a way similar to the present paper. 

§5. Summary 

We have discussed that characteristic aspects of fer- 
romagnetic superconductor UGe2 can be understood in 
a unified way on a single assumption that the supercon- 
ductivity is mediated by the coupled CDW and SDW 
fluctuations originating from the CDW ordering of the 
majority-spin band at T = T^- The growth of extra mag- 
netization below r < Tx is explained by a mode-coupling 
effect of the coupled CDW and SDW ordering caused by 
the imperfect nesting of the majority-spin band. The 
anomalously large entropy around T = Tx(P = 0) is un- 
derstood by the Kohn effect of softened optical phonons 
due to the onset of the CDW ordering of majority-spin 
band. 

We have constructed the general framework of the 
pairing mechanism mediated by the coupled CDW and 
SDW fluctuations. The pairing interaction in the "p"- 
wave channel among quasiparticles in the majority-spin 
band is shown to induce the ferromagnetic fluctuations 
enhanced around P = due to the coupled CDW 
and SDW fluctuations. This mechanism also provides a 
natural explanation for the fact that the superconduct- 
ing phase appears inside the ferromagnetic phase in the 
temperature-pressure (T-P) phase diagram. The anoma- 
lous temperature dependence of the upper critical field 
Hc2{T) observed at P — 11.4 and P = 13.5 kbars is un- 
derstood by the strong-coupling formalism and the ex- 
perimental fact that Tx is an increasing function of H. 

Of course, there still remains many aspects to be clar- 
ified both experimentally and theoretically. A direct 
measurement of coupled CDW and SDW ordering at 
r < Tx is indispensable. In particular, the calculation 



of x(q) based on the band-structure calculation and the 
search of CDW ordering by elastic neutron and X-ray 
scattering measurement is desired especially at around 
P — % kbar where the decrease of p(T) at T = Tx is 
much more prominent compared to the ambient pres- 
sure. Soft modes of phonons toward T = Tx should be 
observed at the wave vector Q where the nesting oc- 
curs. Origin of a first-order like ferromagnetic transition 
at T = Tc has not yet been clarified. The spontaneous 
magnetic moment, in the limit T ^ 0, also exhibits a 
first-order like transition at P = Px This is also to 
be clarified. The identification of the type of the su- 
perconducting order parameter is also important task 
left in the future. To identify the order parameter the 
group-theoretical classification is usually done. However, 
it is necessary to clarify the relevant effective interaction 
which works quasiparticles to identify the most stable 
order parameter among the possible candidates.^' ^-^ We 
suppose that the superconducting gap opens at the parts 
of the Fermi surface with crystalline inversion symmetry, 
which is different from those with the CDW gap to gain 
both the condensation energies of CDW and supercon- 
ducting states. The steep decrease of the superconduct- 
ing temperature, Tgc for the P < Px region in the T-P 
phase diagram may be due to the partial collapse of the 
Fermi surface by the imperfect nesting. 

Quite recently, it has been discovered that ferromag- 
netic metal URhGe, at T < 10 K, sets in superconduct- 
ing state at T ~ 0.35 K under the ambient pressure. 
This compound seems to share common aspects with 
UGe2 , such as the crystal structure with zigzag chain of 
U atoms, existence of a small hump of C/T at T ^ Tc/3 
as seen in Fig. 6, enhanced ratio of 10 times larger 

than the so-called Kadowaki- Woods ratio, suggesting 
that the system is located near the quantum critical point 
in a kind of another, and of course the coexistence of su- 
perconductivity and ferromagnetism. The experimental 
evidence of the CDW has not been reported in URhGe 
at present. However, we expect that the small hump of 
C/T at T - Tc/3 may be the signature of the CDW 
transition, since the zig-zag-chain structure of U atoms, 
which is a common feature of a-Uranium, UGe2 and 
URhGe2, provides the low dimensionality of the Fermi 
surface which enlarges the nesting instability. Further 
studies of URhGe as well as UGe2 are expected to extend 
the frontier of a class of ferromagnetic superconductor. 

Acknowledgments 

This work started from conversations with members 
of Flouquet group at DRFMC of CE A/Grenoble when 
one of us (K.M.) was staying there. We have much ben- 
efited from discussions with A. Huxley, I. Shcikin, D. 
Braithwaite, and J. Flouquet. Unpublished experimen- 
tal data they showed us prior to publication was essential 
for us to develop a theoretical idea. We acknowledge N. 
Tateiwa, T. C. Kobayashi, and Y. Onuki for allowing us 
to use their experimental data of specific heat as well 
as fruitful discussions. We also acknowledge H. Yam- 
agami for providing his result of band structure calcu- 
lations prior to publication as well as valuable discus- 
sions. K.M. has benefited from conversations with G. 



12 



Shinji Watanabe and Kaaumasa Miyake 



Lonzarich, V. P. Mineev, and C. Pfleiderer on the pairing 

interaction, which led us to clarify the presentations in 
§4. S.W. thanks M. Sigrist for useful discussion. Thanks 
are also due to D. Aoki, K. Hori, G. Oomi, N. K. Sato, 
R. Settai and H. Kohno for useful discussions. This work 
was supported by the Grant-in-Aid for COE Research 
(NO.10CE2004) from Monbu-Kagaku-sho, and in part by 
a Grant-in-Aid for Scientific Research for Specified Area 
(No. 12046246) by the Japan Society for Promotion of 
Science. 



quet, J.-P. Brison, E. Lhotel, C. Paulsen: Nature 413 (2001) 
613. 

30) D. Aoki: private communications. 



S. S. Saxena, P. Agarwal, K. Ahilan, F. M. Grosche, R. K. 
W. Hasselwimmer, M. J. Steiner, E. Pugh, I. R. Walker, S. R. 
Julian, P. Moiitlioux, G. G. Lonzarich, A. Huxley, I. Sheikin, 
D. Braithwaitc and J. FLouquet: Nature 406 (2000) 587. 
A. Huxley, I. Sheikin, E. Ressouche, N. Kernavanois, D. 
Braithwaite, R. Calemezuk and J. Flouquet: Phys. Rev. B 
63 (2001) 144519. 

A. M. Clogston: Phys. Rev. Lett. 9 (1962) 266. 

K. Machida and T. Ohmi: Phys. Rev. Lett. 86 (2001) 850. 

L A. Fomin: JETP Lett. 74 (2001) 111. 

N. Tateiwa, T. C. Kobayashi, K. Hanazono, K. Amaya, Y. 
Haga, R. Settai and Y. Onuki: J. Phys. C. 13 (2001) L17. 
G. Oomi, K. Nishimura, T. Kagayama, S. W. Yun and Y. 
Onunki: Physica B 186-188 (1993) 758. 
A. Huxley, private eommunieations. 

N. Tateiwa, K. Hanazono, T. C. Kobayashi, K. Amaya, T. 
Inoue, K. Kindo, Y. Koike, N. Metoki, Y. Haga, R. Settai 
and Y. Onuki: J. Phys. Soc. Jpn. 70 (2001) 2876. 

G. Motoyama, S. Nakamura, H. Kadoya, T. Nishioka and N. 
K. Sato: Phys. Rev. B, in press. 

L Sheikin, private communications. 

L Sheikin, A. Huxley, D. Braithwaite, J.-P. Brison, S. Watan- 
abe, K. Miyake and J. Flouquet: Phys. Rev. B 64 220503 

(2001). 

H. Yamagami and A. Hasegawa; Physica B 186-8 (1993) 182; 
S. Tejima, H. Yamagammi and N. Hamada: preprint. May 
2000. 

A. B. Shick and W. E. Pickett: Phys. Rev. Lett. 86 (2001) 
300. 

H. Yamagami, private communications. 

J. C. Marmeggi, R. Currat, A. Bouvet and G. H. Lander: 
Physica B 263-4 (1999) 624. 

M. Naito and S. Tanaka; J. Phys. Soc. Jpn. 51 (1982) 219. 
R. Settai, M. Nakashima, S. Araki, Y. Haga, T. C. Kobayashi, 
N. Tateiwa, H. Yamagami and Y. Onuki: J. Phys.: Condens. 
Matter 14 (2002) L29. 

N. Suzuki and K. Motizuki: in Structural Phase Transition in 
Layered Transition- Metal Compounds, edited by K. Motizuki 
(Kluwer Academic Publishers, Dordrecht, 1986), P. 135. 
T. Matsuura, K. Miyake, H. Jichu, Y. Kuroda and Y. Na- 

gaoka: J. Magn. Magn. Matter. 52 (1985) 239. 

M. T. Bcal-Monod, C. Bourbonnais and V. J. Emery; Phys. 

Rev. B 34 (1986) 7716. 

K. Miyake, S. Sehniitt-Rink and C. M. Varma: Phys. Rev. B 
34 (1986) 6554. 

D. J. Scalapino, E. Loh, Jr. and J. E. Hirsch: Phys. Rev. B 
34 (1986) 8190. 

A. J. Leggett: Phys. Rev. 140 (1965) A1869. 

H. Maebashi and K. Miyake: Physica B 281&:282 (2000) 526; 

K. Miyake and H. Maebashi: J. Phys. Chem. Solids 62 (2001) 

53. 

E. Felfand and N. R. Werthanier: Phys. Rev. Lett. 13 (1964) 
686; Phys. Rev. B 147 (1966) 288. 

M. Schossmann and E. Schachinger: Phys. Rev. B 33 (1986) 
6123. 

J. Flouquet, A. Huxley, L Sheikin, N. Kernavanois, D. Braith- 
waite, E. Ressouche, S. Watanabe and K. Miyake: J. Phys. 
Soc. Jpn. 70 (2001) Suppl. 14. 

D. Aoki, A. Huxley, E. Ressouche, D. Braithwaite, J. Flou- 



